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Abstract 

We consider mechanisms of directed transport in a ratchet model comprising, besides the external freedom 
where transport occurs, a chemical freedom that replaces the familiar external driving by an autonomous dy- 
namics providing energy input, and an internal freedom representing a functional mode of a motor molecule. 
The dependence of the current on various parameters is studied in numerical simulations of our model. In 
particular, we point out the role of the internal freedom as a buffer between energy input and output of 
mechanical work that allows a temporary storage of injected energy and can contribute to the efficiency of 
current generation. 
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1. Introduction 

Research on directed transport at the micro- 
scopic level has benefitted enormously from the in- 
teraction and mutual stimulation of applied studies 
close to experimentally accessible phenomena and 
theoretical analyses, the paradigm case being the 
concept of Brownian motors [TJ[2J|31[3] as it emerged 
from the discovery of motor molecules generating 
directed transport in the living cell [5] . Notwith- 
standing, the distance between the two endeavours 
remains large. An example is the degree of struc- 
tural complexity, reflected in the number of degrees 
of freedom involved, of the models considered on ei- 
ther side: While realistic molecular-dynamics sim- 
ulations of biomolcculcs typically include the full 
set of hundreds or thousands of nuclear freedoms, 
plus possibly some of the surrounding solvent (for a 
review see, e.g., [7]), the key notion of ratchets has 
been developed considering mere point particles in 
one-dimensional periodic potentials [H ■ 

Efforts are under way on both sides to bring these 
ends closer to one another. On the one hand, in 
molecular biophysics, the concept of functional de- 
gree of freedom [TD1[TT] has been conceived for large- 
scale conformational changes on long timescales 
that are immediately responsible, e.g., for the cat- 
alytic activity of a protein, and can be identified 
on basis of objective criteria such as normal-mode 
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analysis jT2] [TJl [H] to establish a hierarchy among 
the modes from slow to fast and from collective 
(global) to microscopic. On the other hand, ele- 
mentary ratchet models are being endowed one by 
one with additional freedoms to study their possible 
role in directed transport |15j . 

The present work is intended to insert another 
stepstone from the side of abstract ratchet mod- 
els towards biophysical realism, by incorporating a 
few crucial features of molecular motors: We de- 
vise a model that includes an internal freedom to 
represent some functionally relevant conformational 
change. The concept of molecular combustion mo- 
tor is incorporated replacing the familiar periodic 
external driving by a chemical freedom as coher- 
ent energy source which maintains the system far 
from thermal equilibrium [16] . The ubiquitous role 
of thermal fluctuations is modeled, as usually, by 
random forces. At the same time, we avoid the 
detailed reconstruction of any specific system such 
as Fi-ATPase [12 [HI [19] or even their reduction 
to "mechanical toy models" [20] which typically 
still are tailored to represent a certain (species of) 
motor molecule and therefore involve some contin- 
gency. This allows us to study the effect of internal 
freedoms on general grounds, applicable to a broad 
class of molecules and even beyond the biophysical 
realm. In this way, we take up themes set by pio- 
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neering works such as Kaneko et al. |21| on molecu- 
lar machines with internal freedoms and of Mateos 
on ratchets based on non-point particles [15]. At 
the same time, focussing on the biological context, 
we leave aside quantum effects which probably are 
of minor relevance for transport in the living cell, 
yet may be crucial in artificial nanodevices |22j . 

The dynamics of our model is simulated in the 
underdamped regime, allowing for complex deter- 
ministic motion, and in the presence of thermal 
noise, analyzed as to transport mechanisms in- 
volving the internal freedom as essential element. 
We find evidence, above all, that it serves as a 
temporary energy storage which partially decou- 
ples directed transport from the energy source and 
thus decorrelates the discrete steps in the external 
freedom from the likewise approximately quantized 
discharges of chemical energy into the system — 
certainly desirable features from the point of view 
of robustness, efficiency, and versatility of a molec- 
ular motor. Preliminary results of this work have 
been published (in Spanish) in Ref. |23j . 

Our model is motivated and constructed in Sec. [2] 
Section [3] provides a survey of its dynamical be- 
haviour in different parameter regimes. Transport 
properties and the role of the internal freedom are 
discussed in Sec. |U We conclude in Sec. [5] with 
some remarks on open ends. 



standard ratchet potential, periodic with pe- 
riod A ox = 2n but breaking invariance under 
parity x cx — > — x cx , which can be considered ei- 
ther as cyclic (angle) or as extended. To model 
the conversion of internal kinetic energy into 
directed transport, it will be coupled to the 
internal freedom, but also to the chemical free- 
dom from which it receives coherent energy. 

ii. Internal freedom Following Kaneko et al. [21 , 
we model the internal freedom xi n like a pendu- 
lum, with a smooth periodic potential, period 
X\ n = 27r, which however lacks the intentional 
asymmetry of the ratchet potential proper. 

iii. Chemical freedom In order to reproduce the 
entropical bias of the "combustion" reaction, 
we impose a non-zero mean gradient on the 
otherwise periodic dependence of the poten- 
tial on the chemical coordinate x c h, i.e., 

^(^exi -^in: -^ch ~l~ ^eh) ^(^cxi -^in: -^ch) E c \ 1: 

E c h denoting the net energy gain of the (hy- 
drolysis etc.) reaction per molecule. To keep 
the behaviour of this coordinate close to an ap- 
proximately periodic energy injection, even in 
presence of the backaction from the external 
freedom, we further choose a large inertia m c h. 



2. Model 

We pretend to construct a minimal model of a 
molecular combustion motor that goes beyond the 
familiar ratchet scheme in two essential points, (i) 
it includes an internal freedom to represent a func- 
tional (slow, conformational) mode of the molecule, 
and (ii) it incorporates the injection of chemical 
energy, e.g., through hydrolysis of ATP, as au- 
tonomous dynamics of a chemical degree of freedom 
corresponding to the reaction coordinate underly- 
ing that process [16] . The formulation of our model 
is inspired in various respects by the example of Fl- 
ATPase, a prototypical rotational molecular motor 
[17l \TE\ \19\ , but does not intend its reconstruction 
in any biophysical detail, i.e., is to be considered an 
impressionist view, at most, of that molecule. 

More precisely, we require the model to comply 
with the following conditions: 

i. External freedom The coordinate x cx is equiv- 
alent to the single freedom of the usual point- 
particle ratchets. It will be subjected to a 




Figure 1: Potential ([TJ as a function of coordinates x ex 
and x c h for the external and the chemical freedoms, resp. 
Greylevel code ranges from black (low) to white (high). 



To implement these requirements, we compose 
the full potential as follows: 

VijEex: ^iri) ^ch) — ^ei(-^ex) ^in) ~t~ ^n(>^in)~t~ 
+Kc( ) + V r c h(«ch)- 



where 



Vin(^in) 

V C h(a;ch) 



1 + e ci sin(x in ) 

1 +£ei 

sin(cc ex ) + Asin(2x cx ) 

Wa 
Ksm(x in ), 

£ cc sin(x c h) sin(x cx - 5), 
~2n 



Typical parameter values used in our simulations 
are summarized in Tab.[TJ 
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The parameter A determines the asymmetry of the 
ratchet potential [24] . The coupling between exter- 
nal and internal and between external and chemical 
freedoms is controlled, respectively, by e e i and e C c- 
We plot the potential ([!]) as a function of exter- 
nal and chemical coordinates in Fig. [I] Diagonal 
channels corresponding to a rigid association of one 
hydrolysis reaction per spatial step in x ex can be 
discerned (cf. Ref. [US]). 

Moreover, to represent the fast internal as well as 
ambient degrees of freedom in a collective manner, 
we include dissipation terms for x ox and x c h corre- 
sponding to velocity-proportional (Ohmic) friction. 
Working on the microscopic level, we complement 
them by random forces representing thermal noise, 
with autocorrelation functions related to dissipa- 
tion rates through the fluctuation-dissipation the- 
orem 25J. All in all, we thus arrive at a set of 
coupled stochastic differential equations [26 , 



1 + e ci sin(xi n ) 



1 +c ci 

cos(xc X ) + 2Acos(2:c GX ) 
X Wa 

+ £ cc sin(a: ch ) cos(a; ox - S) + £(*), 

TOiniin = KcOs(x in ) 

£ ci cos(x in ) sin(z cx ) + Asin(2a; cx ) 

l+£ci If A 

AV 

m ch x ch = - m ch j ch x ch + 



2tt 

; cos(a; c h) sin(a:e X - 5). 



(6) 



The random force £(t) is defined by = and 

(C(i)C(O)} = 2m cx -/ cx kTS(t) at temperature T. The 
parameters 7 0X and 7 c h are the friction coefficients 
for the external and the chemical freedoms, respec- 
tively. 

The numerical results presented in the following 
have been obtained by solving Eqs. ([6| by Conven- 
tional Brownian Dynamics, see Refs. [56J HE]- 



Parameter 


™ex 


m in 




K 




^ei 


Value 


1 


0.1 


100 


0.1 


1.0 


1.0 


Parameter 


A 


kT 


AV 


5 


7ex 


7ch 


Value 


0.25 


1 


100 


4.53 


0.1 


1 



Table 1 : Default values for the parameters of our model used 
wherever not indicated otherwise. 



3. Dynamics and phase-space structure 

The principal parameters determining the be- 
haviour of our model are A (asymmetry of the 
ratchet potential), e e i (coupling external to internal 
freedom), e GC (coupling external to chemical free- 
dom), 7 ox , 7 c h (friction) and temperature T. This 
quite high-dimensional parameter space is delim- 
ited by the following asymptotics: 

i. Symmetric external potential For A — > 0, 
the simultaneous presence of symmetry-related 
counterpropagating trajectory pairs identically 
cancels transport. 

ii. Internal freedom decoupled In the case £ e i = 0, 
the system reduces to a point-particle ratchet, 
broadly documented in the literature in all 
its facets (overdamped, underdamped with 
nontrivial dynamics beyond mere relaxation, 
Hamiltonian), see [9]. 

iii. Chemical freedom decoupled As we are deal- 
ing with an autonomous system without exter- 
nal driving, decoupling the chemical freedom 
deprives the ratchet of its energy source, so 
that transport dies out. Even in the Hamil- 
tonian case, a driving is necessary for directed 
transport to occur, since there is then no other 
means to break time-reversal invariance (apart 
from magnetic fields which however play no 
role in a biochemical context). 

iv. Hamiltonian dynamics In the limit 7cx,7ch 

of vanishing friction, the system becomes 
a Hamiltonian ratchet with internal freedom. 
However, in the absence of an asymmetric ex- 
ternal driving or a magnetic field, the system is 
then time-reversal invariant and directed cur- 
rents are excluded. 



In the following, we consider in more detail a cru- 
cial aspect of the dynamics, the structure of the at- 
tractor and its connectivity along the lattice as a 
function of the "driving force", i.e., the coupling to 
the chemical freedom, as well as friction and noise 
strength. 





(a) 








Figure 2: Attractor of the system j6j| projected onto the 
plane (p ox ,^cx), for e cc = 1.0 (a) and 7.1 (b). The dotted 
line indicates the surface of section underlying Fig. [7] below. 



Figure [2] shows the attractor, as a converged 
trajectory projected onto the (p ex , x ex )-plaxie, for 
£ cc = 1.0 (panel a) and e cc = 7.1 (b). While for 
low coupling we discern a limit cycle, restricted to 
a unit cell of the potential, for strong coupling it 
has given way to a strange attractor which connects 
across the unit cell boundaries along the entire lat- 
tice. 

These two crucial "crises", though, the emer- 
gence of a strange attractor and its change of topol- 
ogy from isolated to conected, are independent of 
one another. This is evidenced in Fig. [3]}, analogous 
to Fig.[2]but at intermediate values of the coupling, 
where we observe a strange attractor still bounded 
within a unit cell. Conversely, Fig. [3}} shows a limit 
cycle extending periodically along the lattice. 

Besides the fundamental role of thermal fluctu- 
ations in directed transport through the "rectifica- 
tion of noise" outside equilibrium, it affects the dy- 




Figure 3: As Fig.[l]but for e cc = 6.0 (a) and 2.5 (b). 




Figure 4: Threshold value of e* c for the transition from lo- 
calized to extended attractors as a function of the bath tem- 
perature T. 



namics in a more common way by softening attrac- 
tor structures and thus contributing to a merger of 
attractors along the lattice. Figure [4] clearly demon- 
strates how noise tends to lower the threshold £* c 
in terms of the coupling to the chemical freedom for 
the transition to transport along the system. 



4. Directed transport 

Wc here measure directed transport in terms of 
the mean velocity of the external freedom, 

I = (±ex) = lim 7 (x cx (t) - x ex (0)) , (7) 

t— > oo t 

taking averages, e.g., over initial conditions within 
the same basin of attraction or over realizations of 
a stochastic force. The generation of non-zero cur- 
rents in ratchet-like systems depends on the follow- 
ing independent general necessary conditions: 

i. The breaking of all binary symmetries involv- 
ing an inversion of momentum, such as in par- 
ticular time reversal and parity [24 . 

ii. The existence of trajectories extending over 
different unit cells of the potential, or, if the 
external freedom is cyclic, corresponding to ro- 
tational as opposed to librational motion. 

Concerning condition (IB, in the absence of an ex- 
ternal driving or a coupling to a magnetic field, it 
is only the friction terms in the equations of motion 
which break time-reversal invariance. It is evident 
from Fig. [5] that there is no transport for vanishing 
friction 7 CX = 0, and likewise for vanishing asym- 
metry parameter ^4 = (not shown). 

0.12 i 1 1 1 1 1 1 
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Figure 5: Mean current ^fy as a function of the friction coef- 
ficient 7 ex for the external degree of freedom, for e ec = 7.08 
(full line) and 5.35 (dashed). 

As to (JTTJ) , in the presence of friction, the exis- 
tence of persistent transporting trajectories is not 
guaranteed, in difference to the case of Hamilto- 
nian ratchets where there is always an asymptotic 
regime of quasi- free motion at high energies |29| . 
As discussed above, the topology of attractors de- 
pends sensitively on couplings, friction parameters, 
noise strengths etc. In the subsequent paragraphs, 



we consider in detail the influence of a these param- 
eters on the generation of directed currents in the 
system. 

4-1- Resonance between chemical and external free- 
doms 




Figure 6: Mean current ^fy as a function of the coupling e ec 
between external and chemical freedom, for different values 
of the coupling between external and internal freedom, £ e i = 
(full line), 0.5 (bold dashed), 1 (light dashed). 
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Figure 7: Attractor of the system j(S| in terms of the surfaces 
of section indicated in Figs. [2] and [3] as a function of the 
coupling £ ec . Upper part: Mean current |Tf| (bold dashed) 
as a function of the same parameter. 

In the absence of an external periodic driving, our 
model is lacking a precise clock that would lead to 
clearcut resonances. The high inertia of the chem- 
ical freedom, however, and the corresponding rel- 
atively regular injection of energy "quanta" also 
defines an approximate time scale that gives rise 
to resonance-like phenomena. This is obvious from 
Fig. [6j where marked peaks of the current occur at 
approximately equidistant values of the external- 
to-chemical Coupling, £ e c,n — n ^cc,l, n = 1,2,.. 



with e CC! i w 1.8 in this case. For this value, there 
is a 1:1 resonance between the two coupled free- 
doms, i.e., in terms of Fig. [IJ the system is moving 
along diagonal valleys in the two-dimensional array 
of potential maxima. 

However, this resonant behaviour does not 
amount to simple oscillatory motion in a; cx . As is 
evident from Fig. [7J where we plot a section of the 
attractor, cf. Fig. [2j vs. e cc , the peaks of the cur- 
rent coincide precisely with intervals in e cc where 
the attractor is fractal. What is more, the sharp 
definition without shoulders of these peaks can be 
related to abrupt changes in the attractor structure 
from limit cycle to chaotic and back. 

4-2. Effects of noise 
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Figure 8: Mean current ^fj as a function of the bath temper- 
ature T. Inset: enlargement of the low-temperature regime 
< T < 1.5. 

Figure [8] shows the current vs. noise strength in 
terms of bath temperature T . As would be ex- 
pected, too strong noise washes out the structure of 
the attractor, in particular its asymmetry, a neces- 
sary condition for directed transport to occur. The 
inset reveals a more interesting phenomenon: For 
the other parameter values in question kept con- 
stant, there is no extended attractor in the absence 
of thermal fluctuations, so that the current vanishes 
at T — 0. As a consequence, there exists a maxi- 
mum, albeit broad, of the current at a bath tem- 
perature T ~ 1, in analogy to stochastic resonance 

4-3. Role of the internal freedom 

A first indication of the influence of the inter- 
nal freedom can be observed in Fig. |6j comparing 
the width of the current peaks for different values 
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Figure 9: Kinetic energy contained in the chemical (upper 
graph), internal (middle) and external (lower graph) free- 
doms as a function of time. The figure shows a particular 
case. Time dependence is slightly smoothed to suppress irrel- 
evant fluctuations. Arrows indicate typical energy-transfer 
events from chemical to internal (dashed, green) and from 
internal to external freedom (full, red). Small-scale oscilla- 
tions in the lower graph correspond to unit transport steps 
in z ox . See text. 



of the external-internal coupling. They systemati- 
cally grow broader with increasing e c ;. This means 
that off resonance, directed transport is supported 
or even enabled in the first place by the presence of 
the internal freedom. 

Further insight into the underlying mechanism is 
provided by Fig. [9j where we plot the kinetic en- 
ergy in the three respective freedoms as a function 
of time. It shows a clear tendency for a tempo- 
rary storage of energy in the internal freedom: The 
system receives energy almost periodically and in 
nearly equal portions from the chemical freedom. It 
tends to accumulate in the internal freedom (green 
arrow) till a threshold is reached, whereupon it is 
discharged (red arrow) almost completely to the ex- 
ternal freedom, giving rise to a burst of transport 
steps. This "toilet-flush" behaviour is well-known 
from neuron firing [31 j . 

Comparing with Fig. [7J we observe moreover that 
the onset of current with increasing e c j also coin- 
cides with a structural change of the attractor from 
limit cycle to strange. To be sure, two autonomous 
degrees of freedom are sufficient to allow for chaotic 
behaviour. Notwithstanding, we can conclude that 
another relevant role of the internal freedom is to 
induce a richer, more stochastic dynamics of the 
system. 



5. Conclusion 



Acknowledgements 



The present work pretends to reduce the gap be- 
tween abstract models of directed transport and de- 
tailed studies of specific motor molecules by con- 
structing a system that endows ratchets with a few 
essentials of combustion motors, an internal degree 
of freedom and a coupling to an autonomous chem- 
ical freedom as energy source instead of a periodic 
external force, but ignoring the specifics of particu- 
lar (classes of) molecules. This allows us to analyze 
in general terms the role of the internal freedom 
in the generation of transport. We find convincing 
evidence that it serves as an interface which effec- 
tively decouples the output of mechanical work in 
the form of transport steps from the energy input 
through hydrolysis or similar reactions. We can ob- 
serve directly how energy is being stored temporar- 
ily in the internal degree of freedom before being 
channeled to the external freedom. In addition, the 
internal freedom induces transitions towards a more 
irregular dynamics of the external freedom, which 
also contributes to the generation of directed cur- 
rents. 

These results direct attention to processes of en- 
ergy transfer and dissipation within the molecule, 
and in particular indicate the functional advantage 
of a cascaded energy redistribution through inter- 
mediate steps |32j . refining the simple picture of 
an external (transport) degree of freedom coupled 
directly to a structureless bath. Even within the 
scope of our three-dimensional model, such ques- 
tions could be studied by changing the configu- 
ration of couplings, e.g., coupling the external to 
the chemical freedom only through the internal 
one instead of internal-external-chemical as in the 
present work. 

Subtler effects of noisy nonlinear dynamics like 
stochastic resonance [30] have not been considered 
in depth. In the context of our model, it is expected 
to occur in a parameter regime of weak noise close 
to a resonance of the chemical with the external 



freedom, cf. Sec. 4.1 lowering the threshold to di- 
rected transport as a function of the coupling to 
the chemical freedom. Quantum effects have not 
been taken into account, either. We expect rele- 
vant modifications of the picture, if any, for the dy- 
namics of the chemical degree of freedom and the 
injection of chemical energy into functional modes 
where the discretization of energy and quantum co- 
herence might play a major role. 
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